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Abstract 

We study the form factors of the low energy anomalous 7r2y and 737T processes in the non- 
local chiral quark model which incorporates the momentum dependence of the dynamical 
quark mass and realizes correctly the chiral symmetries. The obtained slope parameter for 
7r27 is in reasonable agreement with the direct experimental results but smaller than the 
ones invoking vector meson dominance. Our result for the 737r form factor interpolates 
between the two extremes of theoretical approaches, with the largest one provided by the 
vector meson dominance and the smallest one by the Schwinger-Dyson approach. But all 
of them are well below the single data point available so far. This situation will hopefully 
be clarified by the experiments at CEBAF and CERN. 
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The 7r°77 (vr27) and 77r + 7r°7r~ (737r) processes are the two simplest chiral anomaly- 
driven processes that involve electromagnetic interactions. A consideration of parity con- 
servation, gauge invariance and Lorentz invariance implies the following structures for 
their amplitudes, 

/J7T27 UPt,a Air2y 

Here denote the outgoing momenta of the two photons with Lorentz indices /x and u, 
and p+,o,- the incoming momenta of the three pions, for the two processes respectively. 
The dynamical information is encoded in the form factors A* 21 and A l3n which are Lorentz 
invariant functions of the relevant momenta. In the low energy and chiral limit, they are 
completely determined by the chiral anomaly as summarized in the Wess-Zumino-Witten 
action to be [|l|] 
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where N c and f n are respectively the number of colors and the pion decay constant. 
Beyond the limit, their dependence on the relevant momenta is a reflection of the detailed 
strong dynamics. Since these processes involve only one or a few pions, they may provide 
an ideal testing ground for models of strong interactions. 

The excellent agreement of A^ 1 with the experimental value extracted from the on- 
shell decay of ir° — > 77 had historically constituted one of the first pieces of firm evidence 
that quarks carry three colors. When one of the photons is off-shell, the form factor can 
be parameterized by a slope parameter in the low energy region, 

A n2 yAf = 1 + a x, (3) 

where x = k 2 /™ 2 describes the virtuality of the off-shell photon with momentum k. 
The slope parameter a has been measured both in the time-like region of k using the 
Dalitz decay 7r° — ► e + e~7 and in the space-like region through the ir° production in e + e~ 
collisions. The direct results from TRIUMF and SINDRUM I in the first category are 
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respectively a = 0.026 ± 0.054 § and a = 0.025 ± 0.014(stat.) ± 0.026(syst.) @. The 
CELLO group actually measured the form factor in the large space-like region and then 
extracted the slope parameter by extrapolation using the vector meson dominance to be 
a = 0.0326 ± 0.0026 These results are consistent with each other within the quoted 
errors. Concerning the 737r process the experimental situation is less clear. There has been 
so far one measurement which seems to favor a larger value of Aq 37T than predicted by 
the chiral anomaly. Fortunately this situation will be much improved by the experiments 
at CEBAF || and CERN [|7| which will measure the form factor A 1 ^ in a wider range of 
kinematics. A more precise value of Aq Sw can then be extracted and the form factor will 
be available to distinguish the theoretical results based on hadronic models. 

The low energy physics of the lowest-lying pseudoscalars may be described by a chiral 
Lagrangian which is a tower of terms in increasing order of energy expansion. The struc- 
tures of terms at each order are completely determined by spontaneously broken chiral 
symmetries while their coefficients are left free. These parameters may be modelled by 
properly incorporating the relevant degrees of freedom in the intermediate-energy region. 
Of special interest in this regard are the quark-based models which may have a close 
connection to the underlying QCD dynamics. As is well-known, one feature of dynamical 
quarks is their running mass in the intermediate-energy region, which should have signif- 
icant effects on low energy physics when the quarks are integrated out. This point has 
been nicely taken into account by Holdom and colloborators in their nonlocal constituent 
quark model 00 [0. Indeed, the coefficients in the 0(p 4 ) chiral Lagrangian for the 
lowest-lying pseudoscalars are expressed in terms of convergent integrals of the quark dy- 
namical mass and their phenomenological values are well reproduced. The model has also 
been successful in modelling the low energy hadronic contributions to the running QED 
coupling at the Z boson pole [pT| , and in understanding the quark-hadron duality [|12[ 



and the electroweak couplings of constituent quarks themselves [|13[ . In this note we shall 
examine the other aspect of dynamical constituent quarks, namely, their implications on 
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the anomalous sector of the pseudoscalars, especially the form factors of the 737T and 7r27 
processes . Since the Ward-Takahashi identities for flavor symmetries in QCD are built 
into the model of Holdom et al. we expect that the form factors so obtained should be 
comparable in quality to the coefficients in the 0(p A ) chiral Lagrangian derived from the 
model. Our results will be compared with those based on other approaches. The main 
feature here is that the effects of dynamical quark mass are included in a simplest possible 
form while at the same time avoiding introducing many free parameters. 

The nonlocal constituent quark model is an effective theory in the intermediate energy 
regime. In this model all physics is assumed to be described by a chiral invariant action 
quadratic in quark fields. The dynamical quark mass E(p) is incorporated into the action 
and its momentum-dependent nature leads to nonlocal interactions among dynamical 
quarks, Goldstone bosons and external gauge fields. Let us outline the action involving 
interactions with the external photon field A^x) relevant to our discussion ||10|| . The 
interested reader should consult Refs. [j8| flO[| for a complete account. 



S = J d 4 x fapD^ -jd 4 xj d 4 y E(x - y)$(x)((x)X(x, y)^y)^{y), (4) 

where ip represents the up and down quark fields with dynamical mass whose Fourier 
transform is the quantity And 

Dfj, = d^-ieQA^ Q = diag(2/3,-l/3), 

i = exp(-m 75 // 7r ), 7r = n a T a , 

X(x,y) = Pexp(-iJ?r M (z) dz»), (5) 

T M = i/2[Z{dp-ieQAp)$ + g{dp-ieQA ll )Z] 

= eQA fl + t/(2f^(7Td^-d^7t) + --- } 

where ir a is the pion field, T a is the isospin matrix with Tr[T a T b ] = 5 a t,/2, and P stands 
for path-ordering. For convenience, we list in the following the relevant vertices appearing 
in our calculation of the 737r and 7r27 amplitudes. The QED vertex between quarks and 
the photon is modified to be 

ieQi^-ip + p')^,^)], (6) 
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where p(p') denotes the incoming (outgoing) momentum of the incoming (outgoing) quark 
line (same below), and 

*fe*) = 3^|W (7) 
We should mention in passing that the appearance of the R term in the QED vertex just 
fits the dynamical quark mass £ appearing in the quark propagator so that the Ward 
identity still holds. The pion interaction with quarks is of a familiar form generalized 
from the constant mass case, 

+ (8) 

The model generally contains nonlinear interactions of pions with quarks and photons 
due to the nonlinearly realized chiral symmetry and nonlocality. But we found that for 
the processes considered here only the following interaction involving two pions and two 
quarks can contribute at one loop level, 

^ ({r",7*}[E(p) + E(p + *i)+E(p + fc!,) + E(pO] (g) 
+ X [T a ,T b m P + fc a ) - S(p+ h) + (h - k 2 ) ■ (p + p')R(p,p')]] 



where the two pions carry the isospin indices a, b and the incoming momenta k\, A; 2 
respectively. The parameter x — 0, 1 corresponds to the two versions P||lQfl of the model. 



Since it makes little numerical difference, we shall henceforth take x — 1> corresponding 
to Ref. [fLOfl . As one may easily figure out, only the x term can contribute to Fig. lb for 
the 737T vertex while the first term in Eq. (9) cannot due to symmetry. 

Let us consider the two processes whose Feynman diagrams are depicted in Fig. 1. 
Since we are interested in the form factors in the low energy region, we expand the 
amplitudes in the external momenta. The leading terms must be the same as predicted 
by the WZW action and thus universal to all models which correctly incorporate the chiral 
anomaly. In other words, they must be independent of the specific form of E(p). This is 
indeed the case. For example the leading term in the 7r27 amplitude is proportional to 
the following integral 

d ( x \ 2 
dx— 



o 



dx \x + E 2 ( 1 /!r) 



which is unity independently of S as long as £ is finite in the Euclidean space. For the 737r 
process the leading term is contributed only by Fig. la, whose integral can be simplified 

as 

{ i • ' . ) ./ 1 J ' " 



x + Y, 2 (y/x) J \x + E 2 (a/x) J 
which is always —1/3 for a finite £ in the Euclidean region. The sub leading terms depend 

explicitly on the integrals of £ which are collected using Mathematica. A* 2 " 1 has been 

parameterized in Eqn. (3). For the 737r process, as will become clear later on, we need 

to expand up to the 0(p 4 ) terms to display the kinematic variation of the form factor. 

Using the Bose symmetry we have 

j^y3n 2 6 

— ^ = 1 + m~ 2 £ b % S t + m; 4 £ Ci Q i: (10) 
A) i=i j=i 

where S and Q are symmetrized Lorentz invariants of the momenta pi,2,3 = P+,o,-, 
Si = p 2 +p|+p 2 , 

S-2=P1-P2+P'2-P3+ Pi ■ Pi, 

Qi = (pi? + {plf + (pi) 2 , 

Qi = pIpI + pIpI + pIpI, / u n 

Q-a = P 2 iPi ■ (P2 + Pa) + p\Pi ■ (P3 + Pi) + P3P3 • (Pi + P2), 

Qi = P?P2 • P3 + P|P3 • Pi + P\Pl ■ P2, 

Q5 = Pi ■ P2P2 • P3 + P2 ■ P3P3 ■ Pi + P3 ■ PlPl • P2, 

Qe = (Pi ■ P2) 2 + (P2 ■ Ps) 2 + (P3 ■ Pi) 2 - 
Note that the explicit factors of m, are introduced for convenience although m~ 2 bi and 

m~ 4 Cj actually do not depend on m n . 

The coefficients a, bi and c, are lengthy integrals involving the dynamical quark mass, 

which is in turn related to /V by the Pagels-Stokar formula reproduced in the model 

with T! = £E. A very simple parameterization for S(p) in the Euclidean space was 
suggested by Holdom et al, which incorporates the correct high energy behavior of the 
dynamical mass up to logarithms, 

S(p) = ( A + !) m3 ( 13 ) 
p 2 + Am 2 ' 
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where m is a typical mass scale of the constituent quark and is related to the parameter A 
through the Pagels-Stokar formula. Fixing f n = 84 MeV in the chiral limit we therefore 
have only one free parameter. Since this simple ansatz is quite successful in reproducing 
phenomenological values of low energy quantities as mentioned previously, it will be used 
in our numerical analysis without further adaption. 

Our results for the coefficients a, bi and Cj are presented in Table 1 as a function of 
the parameter A in the same range of values as used previously, where the mass scale m 
is of order 300 MeV. Let us first discuss the slope parameter for the 7r27 process. We 
get a stable result of a = 0.02 for the range of A in the table. This is in reasonable 
consistency with direct results from the Dalitz decays, but smaller than the one extracted 
from the large space-like region by extrapolation using vector meson dominance. The slope 
parameter has been studied in other approaches. The free quark loop |14| with a constant 

2 

constituent mass m predicts a = j^2, which is about 0.014 for m = 330 MeV. In the 
phenomenological approach of vector meson dominance the momentum dependence of the 

2 

amplitude derives from the lowest-lying vector resonances and thus a = ~ 0.03. Chiral 
perturbation theory is appropriate for dealing with low energy pion-photon interactions, 
but it is afflicted in the current case by the unknown counter-term parameters appearing 
in the 0(p & ) anomalous chiral Lagrangian. Assuming they are again saturated by vector 
mesons with a mean mass of my = (9m 2 p + + 2m|)/12, the sum of loop and counter- 
term contributions gives a = 0.032 JT5| . It is clear that our result is larger than the 
one in the constant quark mass model but smaller than the ones (both theoretical and 
experimental) invoking vector meson dominance. 

For the 737r process one has to examine the kinematic variation of the form factor to 
extract information on the coefficients bi and c$. In all of the three experiments available 
or approved, the photon and two of the pions, which we assume to be the first and second 
ones without loss of generality, are on-shell, (pi + P2 + P3) 2 = 0, p\ = p\ = m%- (We 
take m n to be the neutral pion mass below and ignore the small isopin breaking in mass.) 
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Table 1: Results of the coefficients a (in units of 10~ 2 ), bi (10~ 2 ) and q (10~ 3 ) as a function 
of the parameter A. The mass m (in units of MeV) is determined by the Pagels-Stokar 
formula. Ignoring the x term in Eqn (9) would change bi and Cj by less than 10%. 



A 


m 


a 


bi 


b 2 


Cl 


c 2 


c 3 


c 4 


c 5 


c 6 


1 


342 


1.97 


-2.29 


-2.32 


-1.02 


-1.83 


-2.23 


-1.96 


-2.52 


-1.58 


2 


317 


1.94 


-2.37 


-2.32 


-0.97 


-1.77 


-1.95 


-1.71 


-1.89 


-1.22 


3 


299 


1.99 


-2.52 


-2.44 


-1.09 


-1.97 


-2.13 


-1.85 


-1.97 


-1.28 


4 


287 


2.05 


-2.67 


-2.59 


-1.27 


-2.26 


-2.45 


-2.12 


-2.25 


-1.47 


5 


277 


2.12 


-2.82 


-2.73 


-1.47 


-2.59 


-2.83 


-2.43 


-2.61 


-1.71 



The experiment at Serpukhov and the one at CERN are of Primakoff type so that the 
third pion is also on-shell, p 2 = m 2 , while the CEBAF experiment is to be done at a low 
momentum transfer of p\ ~ — m\. Defining the Mandelstam variables s = (pi + P2) 2 , 
t = G°2 + P3) 2 and u = (p3 + Pi) 2 , the form factor is a function of s and t with other 
kinematic variables completely fixed. It is then clear that there is no s or t dependence 
in the 0(p 2 ) terms of A l3n and this is the reason why we expand up to 0(p 4 ). 

We plot in Fig. 2 our numerical results of the form factor A^/Aq 3 ^ at A = 1 as a 
function of s with fixed t = — m 2 , for the Primakoff case (panel (a)) and the CEBAF case 
(panel (b)) respectively. Also shown are the results of other approaches, including the 
free quark loop with a constant constituent mass [IB|, the Schwinger-Dyson approach in 
the generalized impulse approximation Jl7| , chiral perturbation theory with vector meson 
saturation ||18|| , vector meson dominance JL9] and its unitarized version [20|. The form 
factors expanded up to second order in s and t in the free quark loop and the Schwinger- 
Dyson approaches can be read off in the original papers. The chiral perturbation result 



augmented with vector meson saturation of counter-terms is [pL8 



Aq 37T 



2m 2 

5 . . Ami 

+ -(s + t + u) + ^ 



1 / s , Tni 

- (s + t + U lB-2 

3 mi 



f(m 2 n ,s) + f(ml,t) + f(ml,u) 



(14) 
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where 



/(m 2 ,g 2 ) 



z + 1 

'1 — x)z In 2, for x < 0, 

z — 1 

^1 — x)^ 2 arctan 2, for < x < 1, 

In Z7T 



1 — 12 



1 



2, for 1 < x, 



(15) 



.r 



1 - 

4m 2 ' 



.r 



which is different from the one quoted for the Primakoff case in Ref. [20[|. The phe 



nomenological approach of vector meson dominance gives |i9 

,2 



which is unitarized to be pO 



D!(q 2 



1 



1 - 



+ 



^ m 2 — s m 2 — t 



+ 



(16) 



+ 



+ 



m. 



(m 2 — s)(m 2 — t)(m 2 — u) 



rn 



96vr 2 / 2 m 2 24vr 2 / 2 



(17) 



Note that the results for chiral perturbation and unitarized vector meson dominance are 
actually shown for |A 737r /y4o 3?r | since the form factor can become complex in these cases. 

It is clear from Fig. 2 that the Schwinger-Dyson approach always gives the lowest 
values of the form factor while the vector meson dominance (especially its unitarized 
version) predicts the largest values and the steepest change in the kinematic region con- 
sidered here. It is interesting to notice that in contrast to the case of the vertex tt2^ the 
chiral perturbation theory predicts a much lower value of the 737r amplitude than the 
vector meson dominance does. Our results interestingly interpolate the two extremes and 
are slightly larger than the one using a constant quark mass of 330 MeV. 

We have studied the form factors of the low energy anomalous 7r27 and 737r processes 
in a simple quark-based model which incorporates the momentum dependence of the 
dynamical quark mass and realizes correctly the chiral symmetries. The obtained slope 
parameter for 7r27 is in reasonable agreement with the direct experimental results from 
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TRIUMF and SINDRUM but smaller than the ones (both theoretical and experimental) 
invoking vector meson dominance. All theoretical predictions for the 73-7T form factor are 
well below the single data point available so far. But there are also significant differences 
among these theoretical results. This situation will hopefully be clarified and distinguished 
by the experiments at CEBAF and CERN. 
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Figure Captions 



Fig. 1 Feynman diagrams for the vertices jStt (a and b) and tt27 (c). The solid, dashed 
and wavy lines stand for the quark, pion and photon fields respectively. 
Fig. 2 The form factor A^ 3n /Af n at A = 1 (solid curve) as a function of s/m^ (m n = 135 
MeV) for the Primakoff case (panel (a)) and the CEBAF case (panel (b)) respectively. 
Also shown are the results of the following approaches: the free quark loop with a con- 
stant constituent quark mass of 330 MeV fllEfl ; the Schwinger-Dyson approach |T7j]; chiral 
perturbation with vector meson saturation [I8|; vector meson dominance ]Tjl and its 



unitarization 20 
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Figure 2 
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